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(T) ' Abstract. We present dust features and masses observed in young supernova 

remnants (SNRs) with Spitzer IRS mapping and staring observations of four 
youngest supernova remnants: SNR 1E102.2-7219 (E0102) in the SMC, Cas A 
and Gil. 2-0. 3 in our Galaxy, and N132D in the LMC. The spectral mapping 
data revealed a number of dust features which include 21 micron-peak dust and 
■ featureless dust in Cas A and 18-micron peak dust in E0102 and N132D. The 

(— | | 18 micron-peak feature is fitted by a mix of MgSiC>3 and solid Si dust grains, 

O j. while the 21-micron peak dust is by a mix of silicates and FeO; we also explore 

I ■ dust fitting using Continuous Distribution of Ellipsoid grain models. We report 

detection of CO fundamental band from Cas A in near-infrared. We review dust 
features observed and identified in other SNRs. The dust emission is spatially 
correlated with the ejecta emission, showing dust is formed in SN ejecta. The 
spectra of E0102 show rich gas lines from ejecta including strong ejecta lines of 
Ne and O, including two [Ne III] lines and two [Ne V] lines which allow us to 
diagnostic density and temperature of the ejecta and measure the ejecta masses. 
E0102 and N132D show weak or lacking Ar, Si, and Fe ejecta, whereas the 
young Galactic SNR Cas A show strong Ar, Si, and S and weak Fe. We discuss 
compositions and masses of dust and association with those of ejecta and finally, 
dust contribution from SNe to early Universe. 
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OS . 1. Introduction 

o 

Meteoritic and astronomical stuides show that presolar, cosmic grains condense 
in the dense, warm stellar winds of evolved stars and in supernova explosions. 
Mantles of the pre-existing dust in molecular clouds are vaporized as the forming 
stars and planetary system heat them. A small fraction of the dust survived 
solar system formation without alteration, protected inside asteroids. The most 
abundant presolar grains are SiC, nanodi amonds, amorphous si licates, forsterite 
and enstatite, and corundum (AI2O3) ([Messenger et al.l 120061 ). Some isotopic 



anomalies of heavy elements in meteorites have been attributed to the dust that 
had condensed deep within expanding supernovae. 

Recent deep sub-mm observations have also shown there to be galaxies 
and QSOs with very large dust masses (>10 8 M Q ) at z> 5 (Bertoldi et al. 
2003). The timescales for low-mass stars to release their dust are too long to 
explain these high redshift systems. In contrast, supernovae produce copious 
amounts of heavy elements and release them on short timescales (Morgan &: 
Edmunds 2003; Dwek et al. 2008). Theoretical modeling of the conditions in 
the supernova ejecta indicate that Type-II SNe are sources of dust formation and 
should be able to pr oduce substantial quantities of dus t , on the order of a solar 
mass per explosion ( Deneault. Clayton. Sz Heger 1120031 ; Todini Ferraral 12001 



Nozawa et ahll2003l h Yet until very recently there existed little observational 



evidence that this actually occurs. 

There is now clear evidence for dust formation in core-collapse supernovae, 
but the quantity of dust formed within SNe ejecta is still a subject of debate. 
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For SN1987A, this includes dust emission, dust absorption and a drop in line 
intensities for the refractory elements that signals that dust is being formed 
(Kozasa et al. 1989). Detections of SiO and CO fundament als with Spitzer ob - 
servations futher provide evidence of dust formation in SNe ( Kotak et al. 2006). 
For the young Galactic SNR, CasA, ISO observations gave evidence for the as- 
sociation of the dust with the ejecta by requiring a mixture of dust grains that 
are not t ypical of the ISM. Submillimeter observa tions of Cas A and Kepler with 
SCUBA (|Dunne et al.ll2003l : iMorgan et al.l 120031 ) suggest the presence of large 
amounts of cold dust (~ 0.3 — 2M ff at 15-2 K) but with some controversy 
related to foreground material. Krause et al. ( 20041 ) showed that much of the 
160/im emission observed with Spitzer is fo r egrou nd material, suggesting there 



is no cold dust in Cas A. [Wilson Sz Batrla I ( 20051 ). however, used CO emission 



towards the remnant to show that up to about a solar mass of dust could still 
be associated with the ejecta rather than with the foreground material. 

We present Spitzer IRS spectral mapping of the entire SNRs of Cas A, 
E0102, and N132D and an IRS staring mode data of Gil. 2-0. 3, and near-infrared 
observations of Cas A for the first overtone CO detection. We identify the dust 
composition associated with the ejecta and discuss dust formation in SNe with 
an estimate of the total mass of freshly formed dust. 



2. Cas A 

Cas A is one of the youngest Galactic SNRs, with an age of 335 yr attributed to 
a SN explosion in AD 1671. The proge nitor of Cas A is b elieved to be a Wolf- 
Rayet star with a mass of 15-30 Mq ( Young et al. 20061 ) We present a short 
summary of Spitzer Infrared Spectrograph (IRS) mapping observations of Cas 
A which were published in Rho et al. (2008), and we added spectral fitting with 
the dust models of continuous distributions of ellipsoidals (CDE). 
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Figure 1. (a: left) The 21 /im dust map of Cas A shows remarkable sim- 
ilarity to the Ar ejecta map, showing that dust is formed in the ejecta. (b: 
top right) Representative set of Spitzer IRS spectra of Cas A with rich ejecta 
lines, The dust feature at 21 ^m often accompanied by a silicate emission 
feature at 9.8 fim with strong Ar lines (red, curve a) and weak-21 /im dust 
with relatively strong Ne lines compared with Ar (green curve, b). Feature- 
less spectra include the continuous rising spectra (c and d). (c: bottom right) 
Featureless dust Spectrum: the continuum can be fit with MgSi03 and Fe 
accompanied with S, Si, and O/Fe lines. 
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Gas lines of Ar, Ne, 0, Si, S and Fe, and dust continua were strong for 
most positions. We identify three distinct ejecta dust populations based on their 
continuum shapes. The dominant dust continuum shape exhibits a strong peak 
at 21 /im. A line-free map of 21 /im-peak dust made from the 19-23 /im range 
closely resembles the [Ar II], [O IV], and [Ne II] ejecta-line maps as shown in Fig. 
[Th., implying that dust is freshly formed in the ejecta. We identify three distinct 
classes of dust associated with the ejecta as shown in Fig. [lb. Spectral fitting 
of the 21/im-peak dust implies the presence of Si02, Mg protosilicates, and FeO 
grains in these regions in Fig. [2k. The silicate composition is responsible for the 
21 /im peak, suggesting that the dust forms around the inner-oxygen and S-Si 
layers and is consistent with Ar being one of the oxygen burning products. We 
also include amorphous MgSiOs (480 K) and SiC>2 (300 K) to account for the 
emission feature around the 9.8 /im. The composition of the low temperature 
(40-90 K) dust component necessary for reproducing 70 /im is rather unclear. 
Either A1 2 3 (80 K) (Model A in Table 1 of Rho et al. 2008) or Fe (100 K) 
(Model B in Table 1) can fit equally well. 

The second dust type exhibits a rising continuum up to 21 /im and then 
flattens thereafter. This "weak 21 /im" dust is likely composed of AI2O3 and 
C grains. The third dust continuum shape is featureless with a gently rising 
spectrum and is likely composed of MgSiO"3 and either AI2O3 or Fe grains (Fig. 
[Tfc) . Using the least massive composition for each of the three dust classes yields 
a total mass of 0.02 M . Using the most-massive composition yields a total 
mass of 0.054 M . The primary uncertainty in the total dust mass stems from 
the selection of the dust composition necessary for fitting the featureless dust 
as well as 70 fim flux. The freshly formed dust mass derived from Cas A may 
be able to explain the lower limit on the dust masses in high redshift galaxies 
when assuming there is no significant dust destruction. 

We also performed spectral fitting of the 21/im-peak dust with CDE models 
for SiO"2 and MgSiO"3 and the spherical model for the rest of dust composition. 
The 21/im peak dust is reasonably well fit by with SiO"2 and AI2O3 in Fig. [2b 
without the composition of Mg proto-silicate and with much less contribution 
of FeO. In addition, previously unidentified 11-12.5/iin feature can be fit with 
the CDE model of SiC, which is one of the important pre-solar grains. Slight 
excess of model at 10/im absorpt i on fea ture can be compensated when we use 
the extinction of Chiar &: Tielen I ( 20071 ) in Fig. [2b- The estimate of dust mass 
using the CDE fittting of Si02 slightly increases the mass, but within the range 
of the previous estimate of Rho et al. (2008). 



3. SNR 1E0102.2-72.3 (E0102) 

We present Spitzer IRS and IRAC observations of the young supernova rem- 
nant E0102 (SNR 1E0102. 2-7219) in the Small Magellanic Cloud. We per- 
formed an IRS staring observation toward the southeastern shell of E0102 (R.A. 
01 h 04 m 04.04 s and Dec. -72°02'00.5", J2000) as a part of our Young SNR Spitzer 
GO program (PI: Rho). The Long Low (LL: 15-40 /im) IRS data were taken on 
2005 August 14 with 6 cycles of 30 sec exposure time; this yields a total exposure 
time of 360 sec for the first and second staring positions. The Short Low (SL: 
5-15 /im) IRS observations were made with 3 cycles of 60 sec exposure time; this 
yields a total exposure time of 360 sec per sky position. We also used archival 
IRS mapping data from the IRS Legacy SMC program, which data were 6 times 
shallower than our GO IRS staring mode data, but still helpful to generate line 
maps. 

The infrared spectra in Figure [3] show both line and continuum. The de- 
tected lines include strong ejecta lines of Ne and O, including two [Ne III] lines 
at 15.5 and 36.0 //m, and two [Ne V] lines at 14.3 and 24.3 /im and weak Si and 
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Figure 2. (a: left) 21 ^m-peak dust spectrum superposed on the dust fit of 
Model A: A Spitzer IRS spectrum towards a bright part of the northern shell 
fitted with dust compositions of Mg proto-silicate, MgSiC-3, SiC>2, FeO, and 
AI2O3. The compositions suggest that the dust forms around inner-oxygen 
and S-Si layers, (b: right) Spectral fitting with continuous distribution of 
ellipsoidal models of SiC>2, MgSiC>3 as well as other composition of AI2O3, 
FeS and FeO, showing CDE model of SiC>2 can reasonably well reproduce the 
21/im-peak dust without the Mg proto-silicate. 



S lines. Unlike the young Galactic SNR Cas A, E0102 lacks emission from Ar 
and Fe, both in the infrared and at other wavelengths. The [Ne II] line at 12.8 
fim is broad, with a velocity dispersion of 2,000-4,500 km s _1 , showing that it 
originates in fast-moving ejecta. 




Figure 3. (a) Spitzer IRS spectrum of E0102. (b) E0102 IRS dust spectrum 
superposed on a dust fit (Model A). The continuum was fitted with dust 
compositions of Si, MgSiC>3 and Carbon. The compositions suggest that the 
dust forms around carbon-burning layers. The data and the total fit are shown 
in blue and thick red lines, respectively. Contributions from each type of dust 
are shown with the dust temperatures given in parentheses. The gas lines are 
also shown in the plot, but they are excluded in spectral fitting 



3.1. Dust Spectral Fitting and Dust Mass of E0102 

The IRS spectra of E0102. 2-7219 have a prominent dust feature peaking at 
18 /im which coincides spatially with the emission from the infrared-emitting 
ejecta. This is clear evidence that dust is forming in the ejecta of the supernova 
remnant E0102. To determine the dust composition and mass, we performed 
spectral fitting to the IRS dust continuum using the deep IRS spectra shown in 
Figured The spectrum shows a 18/iin-peak in the continuum which we attribute 
to dust. First, we estimated the contribution of the synchrotron emission to the 
infrared spect rum using the rad io fluxes (0.65 Jy at 408 MHz) and spectral 
index (-0.70) (jAmv fe Ball 1993). The expected synchrotron fluxes are (1.4 - 
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5.5) xlO -4 MJy sr" 1 between 5-40 //m, which is a very small contribution to the 
infrared continuum (<1%, see FigSJ). 

We performed spectral fitting to the IRS continua with the Planck function 
B V (T) multiplied by the absorption efficiency (Q a bs) for various dust composi- 
tions and the detailed dust fitting technique and mass estimation method is de- 
scribed in iRho et al.l ( 20081 ). The dust compositions, temperatures, and masses 
obtained from the best fits are summarized in Table [TJ 



Table 1. Properties of Freshly Formed Dust in E0102 

Model (Ax^) Dust Compositions (temperature [K], mass [Mq]) 



Mass (M Q ) 



0.015 ( a 5.12xl0- il ) 
0.007 ( a 2.53xl0~ 3 ) 

0.014 ( a 4.80xlQ- 3 ) 



A (2.41) MgSi0 3 (116, 3.80E-5), Si (150, 6.20E-5), C (55, 5.00E-3) 

B (1.98) MgSi0 3 (117, 3.21E-5), Si (114, 4.4E-4), Al 2 :i (=55, 2.06E-3) 

C (2.16) MgSi0 3 (124, 2.34E-5), Si (149, 6.3E-5), 

A1 2 Q 3 (60, 3.33eE-4), C (55, 4.350E-3), Mg 2 Si0 4 (100, 1.10E-5) 



The 18/mi-peak dust feature is fitted by a mix of MgSi03 and Si dust grains, 
while the rest of the continuum requires either carbon (Model A; Figure H|) or 
AI2O3 grains (Model B). The temperature of AI2O3 cannot be constrained by 
our spectra, which lack long-wavelength data, so we set the same temperature 
as the carbon in Model A. We favor carbon or AI2O3 over solid Fe dust, because 
we expect carbon or AI2O3 dust to be present where the Ne and O ejecta lines 
are dominant: N e, Mg, and Al are all carbon-burning nucleosynthesis products 
riWooslevl[l1)95T ). MgSi0 3 and A1 2 3 are condensed in the ejecta of Ne, Mg, 



and Al, which are also carbon-burning products; the O and Al are found in the 
outer layers of ejecta, which is where Si and carbon dust grains condense. The 
dust and ejecta compositio ns seen in E010 2 are similar to those of the weak 
21 /mi dust layer in Cas A ( Rho et all 120081 ) . To explore the effect of including 
many minerals, in Model C (Figure U]) we performed fit including MgSiC>3, Si, 
AI2O3, carbon and Mg2Si04 (the most commonly expected dust compositions 
from Ne and O ejecta; see e.g.TOl, N03). Note that presolar AI2O3 grains from 
meteorites have been inferred to b e among most abund ant isotopically-enriched 
materials ejected by Type II SNe ( Clayton et al.l [20041 ). The mass is relatively 
well constrained by the observations (within 20-40%). 

The estimated dust masses within the IRS observed slit are 5.0xl0 -3 and 
5.5xl0 -3 M for Models A and D, respectively. To estimate the total dust mass, 
we correct a factor of 2.9, the observed (within the high-sensitivity spectral slit) 
mass for the fraction of the SNR covered by the IRS slit using a continuum 
map. The total dust masses for the entire SNR are then 0.015, and 0.014 M© 
for Models A, and D, respectively. Our estimation exceeds previous estimates 
(8xl0 -4 M ; Stanimirovic et al. 2005) by more than an order of magnitude 
because we have accounted for dust compositions. 
3.2. Ne Line Diagnostics and Shock Models 

The ions Ne III and Ne V ions offer pairs of lines in the IR that are suitable 
diagnostics for the density and temperature in emitting regions. The measured 
line flux ratios are [Ne V] A 14.3/24.3/xm = 1.76±0.11, and [Ne III] A 15.6/36/xm 
= 54.3±27.4. To constrain temperatures and densities, we calculate the line 
intensities and ratios of [Ne V] and [Ne III]. We solve the level population equa- 
tions for collisional excitation as a matrix using 5 energy l evels for Ne V and 
3 leve ls for Ne III. The input atomic data were taken from iGriffin &: BadneTT 



(2000). The line diagnostics are shown in Figured The electron density and 
temperature jointly obtained from the ratios of [Ne V] and [Ne III] are 7xl0 4 
cm -3 and ~610 K, when assuming that both Ne III and Ne V gas come from 
the same gas. Such a hi gh density seems to be chara cteristic of SN ejecta as 
observed in optical lines; Chevalier Kirshner I ( 19781 ). for example, suggested 



that the optical [S II] line in Cas A originated from ejecta with a density of 
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Figure 4. (a: left) The dust continuum superposed on a combination of dust 
models of Si, MgSiC>3, AI2O3, Mg 2 SiC>4, and carbon, (b: right) Line diagnos- 
tic contours of the ratios of [Ne III] 15/35^m (red) and [Ne V] 14.3/24.3/xm 
(blue). The observed ratios are marked with thick solid lines for the best 
value of the ratios [Ne III] = 54 and [Ne V] = 1.76. Note that the density 
is an electron density. The larger shaded region shows the range of tempera- 
tures and densities allowed by errors for the [Ne II] ratio, whereas the allowed 
physcial conditions of [Ne V] lines are marked with thick lines. 



10 5 cm -3 . iBlair et~aL (2000), however, find considerably lower densities when 
matching optical lines observed in bright optically emitting knots in E0102. 

We run radiative shock model with a velocity r ange of 50 - 500 km s~ 3 , 
and use basic parameters quoted in Blair et al. (2000) and Hughes et al. ( 19981 ) 
including the abundances. To calculate the hot, non-radiative region of the 
shock, we calculate the steady shock profile, integrating the energy equation and 
including the magnetic field, radiative cooling and non-equilibrium ionization. 
The cooled photoionized region of the shock requires a code that can calculate 
the thermal balance at the low te mperatures an d high densities of that gas and 



for this purpose we use Cloudy (|Ferlandl[200rJ l. We use the EUV/X-ray flux 



generated in the hot post-shock region as the input ionizing flux and use the 
results from Cloudy on the Ne ionization to explore the parameters needed to 
match the observed [Ne II], [Ne III], and [Ne V] lines. The shock model indicates 
that the [Ne II] and [Ne III] lines are mainly from photoinization zone with a cold 
temperat ure of 400 - 1000 K , while the [Ne V] comes from hotter density region. 
Following IBlair et al.1 (|200Ch we calculate a series of shocks with the same ram 
pressure but shock speeds ranging fro m 50 - 500 km s -1 such that nQV^ hock = 



-1\2 



cm )(100 km s" 
field strength of Bq 

114 



Also following IBlair et al.l (|200d ) we assume a magnetic 
4no(cm -3 ) 1 / 2 /xG. The derived [Ne II] column densities of 
(2-40) xlO 14 cm -2 are a factor of a few larger than those of [Ne III] and [Ne I], 
and the total infrared Ne ejecta (from [Ne I] to [Ne V]) are a few 10~ 3 - 10~ 2 
Mq. The uncertainty in the Ne mass is due to uncertainties in the pressure and 
magnetic field in ejecta clumps. 



4. N132D 

N132D is the youngest SNR known in the Large Magellanic Clouds (LMC) 
with a diameter of 100" (physical size of 25 pc) and an age of ~ 3200 yr. 
N132D shows two types of emission, ejecta and ISM emission. Optical ob- 
servations identified O-rich ejecta (like Cas A) emission, an inner ring mor- 
phology, which shows elevatered abundances and red and blue-shifted oxygen 
eject a with a velocity as high as 4400 km s _1 (Mo rse. Winkler, Kirshner 
1995). Ejecta are oxygen-rich like in Cas A. The progenitor is suggested to 
be Type lb, a Wolf-Rayet (W/O) star with a mass between 30 and 35 M 
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Figure 5. (a: left) N132D Spitzer images: MIPS 24/Ltm, [Ne II] (red), [O 
IV] (green), and color composite (red: [Si II]), showing that Ne, O, and Si 
maps are dominated by ejecta, and coincides with the optical O-rich ejecta 
(M orse. Winkler, fc Kirshner 1 11995). The contours show the location of op- 
tically emitting O-rich ejecta superposed on [O IV] map. (b: top right) IRS 
spectrum of N132D (red) compared with that of E0102 (blue), showing re- 
markably similar ejecta composition and dust continuum shape, (c: bottom 
right) The dust continuum superposed on the best fit of dust model, a com- 
bination of MgSiC>3, solid Si and carbon. 



(iBlair et alJl200nh. Recently. X -ray observations identified highly ionized oxygen 
ejecta ( Borkowski. Hendrick. &: Reynolds 20071 ). The bright shell is dominate d 
by shocked ISM, particulary bright in the southeastern shell (jBlair et al.l 12000). 
Our previous Spitzer observations revealed PAH emission from the dominant 
component of 15-20 PAH hum p which is interpreted as PA Hs with a relatively 
large number of carbon atoms ( Tappe. Rho &; Reach I [20061 ). 

We made follow-up observations (PI: Tappe, PID: 30372) with Spitzer for 
IRS-LL mapping covering the entire SNR, and IRS staring of SL observations 
toward 4 positions. Staring mode was 12 cycles of 60 sec exposure for SL. IRS- 
LL mapping is 2 cycle of 30 s exposure for 10" stepsize, resulting in 4 minutes 
integration per position. The observations took place on 2007 Jun 23. 

We present detection of infrared ejecta in N132D as shown in Fig. [5^i. 
The line maps of [Ne III] at 15.5^m and [O IV] 26 fim, and [Si II] 34.8^m 
maps show strong emission at the central ring, which coincides with the known 
optical and X-ray ejecta. Fig. \5jp shows comparison between N132D and E0102 
ejecta spectra; showing remarkably similar lines of Ne and O and the continuum 
shapes. The lines at the ejecta positions are a few ten times stronger than those 
in the shocked ISM position, and N132D shows a higher ratio of [Ne II] /[Ne III] 
than that of E0102, indicating that Ne is less ionized. We performed spectral 
fitting to the IRS dust continuum as for E0102 (see Fig. [3]). The spectrum in 
Fig. \5jp was fit best with the composition of MgSi03 and solid Si and C, and 
the derived dust mass for the ejecta emitting within one IRS slit is 8(=L3)xlO~ 3 
Mq. We can also fit the spectrum with a combination of MgSiOs, Si, carbon 
and/or AI2O3, but the mass could not be well constrained. This will be a lower 
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limit of freshly formed mass, since the rest of the ejecta and fresh dust has been 
mixed with shocked ISM material. 



5. Gil. 2-0.3 



Gil . 2-0.3 is the remnant of the historic supernova of AD 386 ([Reynolds et al 



1994) . An X-ray e mitting 65 ms pulsar is located at the center of the remnant 



(jKaspi et ai1l200ll ) with an age of 1600 yr. The pulsar wind nebula shows elon 



gated, hard X-ray structure s. The distance is 4.4 kpc, based on H I absorption 
( Radhakrishnan et al. 19721 ). Radio observations show bright radio emission, a 



circula r shell of mean di ameter of 4.4', 22 Jy at 1 GHz with an spectral index 



of 0.6 (|Green et al.lll988l ). 

We observed IRS spectra in staring mode toward two positions (18 h ll m 34.76 s , 
Dec. -19°26'30.0"; 18 h ll m 32.32 s , -19°27'10.5"). The observations took place 
on 2007 October 4, and the Long Low (LL: 15-40 /im) IRS data were taken with 
3 cycles of 30 sec exposure time; this yields a total exposure time of 180 sec 
for the first and second staring positions. The Short Low (SL: 5-15 [im) IRS 
observations were made with 8 cycles of 14 sec exposure time; this yields a total 
exposure time of 224 sec per sky position. Reference positions were taken with 
IRS mapping toward nearby positions of ~30' away. We also conducted near- 
infrared observations using Palomar telescope using narrow-filters of Fe and H2 
images. The observations were conducted on 2007 October 22 and October 18. 

Fig. [6]shows the IRS spectrum of Gil. 2-0. 3 towards Fe-rich knots, showing 
strong [Fe II] and Ne lines, and weak Ar, Ni, Si and S lines. Diagnostics of 
the [Fe II] lines use excitation rate equa tions which are presente d in Rho et 
al. (2001) and we updated atomic data ( Ramsbottom et aill2007l ). Figure [6b 



shows contours of line ratios of 17.9/5.35/um and 17.9/26/um; the former ratio 
is mainly sensitive to density and the latter ratio changes depending on both 
density and temperature. The two ratios imply that Fe emission of Gil. 2- 
0.3 has a temperature of ~80,000 K, and an electron density of 6000 cm~ 3 ; 
such as high density is a characteristics of ejecta. In compariso n, the [Fe II] 



line i mplied electron density from the SNR Kes 17 is ~100 cm 3 ([Hewitt et al 
2009). Fe is a rare element in the circumstellar medium, but it is rich element of 
nucleosynthetic yields. The continuum of Gil. 2 shows smooth featureless shape 
indicates possible Fe dust. Detection of Ni lines also implies that they are ejecta 
from nucleosynthetic yields. 

6. Carbon Monoxide in the Cassiopeia A Supernova Remnant 

We report the likely detection of near-infrared 2.29 /iin first overtone Carbon 
Monoxide (CO) emission from the young supernova remnant Cassiopeia A (Cas 
A). The continuum-subtracted CO filter map reveals CO knots within the ejecta- 
rich reverse shock. We compare the first overtone CO emission with that found 
in the well-studied supernova, SN 1987A and find ~30 times less CO in Cas 
A. The presence of CO suggests that molecule mixing is small in the SN ejecta 
and that astrochemical processes and molecule formation may continue at least 
~ 300 years after the in itial explosion. This section is a summary of a Journal 
paper ([Rho et alJl2009al ). 



Introduction and Observations Supernova e (SNe) are suggested to be the 
first molecular factories in the early Universe ( Cherchneff Sz Lilly 2008). Un- 



derstanding the chemistry within the ejecta, and the abundance and mixing of 
CO and SiO molecules, is particularly important in studying dust formation 
processes since this controls how much carbon is available to form amorphous 
C grains and how much oxygen available beyond that bound in CO and SiO 
molecules to form MgSiOa and Mg2Si04. 
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Figure 6. (a: left) The IRS spectrum of Gil. 2-0. 3 (red) compared with 
that of Kes 17 (green: from Hewitt et al. 2009), a typical ISM dominated 
SNR, showing different ratios of 5.3, 17.9 and 26 /im [Fe II] lines, (b: right) 
Contour plots of of [Fe II] diagnostic line ratios using 5.3, 17.9 and 26 /jm 
lines. The observed ratios of Gil. 2-0. 3 (filled dot) and Kes 17 (diamond) are 
marked, showing that Fe emission from Gil. 2-0. 3 has a higher density and 
temperature. 



We observed Cas A on 2006, September 3 and 4 with the Wide field InfraRed 
Camera (WIRC) on the Hale 200 inch (5 m) telescope at Mount Palomar. We 
took narrow-band images using filters centered on CO(2-0) at 2.294 /xm, K- 
band and H-band continuum at 2.270 and 1.570 /im, and P/3 at 1.182 fim (see 
Figure [7|). The exposure time of CO, K-cont and P/3 is a has a total on-source 
integration time of 8100 sec for CO, 7290 sec for K-cont and 810 sec for P/3 
image. 

Results and Discussion Fig. [7^, shows a color composite image of Cas A, 
combining the P/3 line (blue), 2.27 ^m continuum (green), and CO (2-0) 2.29 p 
line (red). The 2.27 /mi continuum consists mostly of non-thermal synchrotron 
radiation and delineates the shell structure (see Rho et al. 2003). Finally, the CO 
(2-0) overtone emission traces the dense knots of gas along the shell and interior 
interface. The molecular emission is evident as red in Figure [7] and more easily 
seen in the continuum-subtracted CO image, revealing the line emission towards 
the northern parts of bright rings and an eastern filament of the Minkowski 
knots in the southeast. The CO emission is highly clumpy and distributed in 
a collection of knots. The locations where CO emission is detected coin cide 
with the ejecta at the reverse shock dHwang et alJ[2003 : iRho et all 120081 ). It 
indicates that the CO gas was formed in the ejecta. We compiled spectral 
energy distributions (SEDs) of near-infrared H, K-continuum, CO bands and 
Spitzer IRAC Bands 1, 2 and 3 (3.6 - 5.8 /im) for bright knots; Fig. [7] shows 
three representative SEDs after correcting for the extinction. The SEDs show 
a sudden jump at 2.294 /im toward the CO detection positions, with no such 
jump toward the NW ejecta position. The SEDs of the CO knots also show a 
larger IRAC Band 2 (4.5 /im) excess than that of the NW ejecta knot. 

We estimated the mass of CO in Cas A by comparing the CO flux density 
in Cas A with that found in the supe rnova SN 1987A. F ig. [7] compares the CO 
emission from Cas A and SN 1987A (jMeikle et al. 19891 ). After accounting for 
the distance differences, the CO (2-0) first overtone flux density of Cas A is a 
factor of ~ 28 smaller than that of SN1987A on day 255. Assuming the same 
gaseous conditions as seen in the CO emitting regions in SN1987A (a tempera- 
ture of 1800-2800 K and a velocity of 2000 km s _1 ), we estimated the CO mass in 
Cas A to be M C o ~ (1.7-8.2) x 10" 6 M . The CO detections in SNe have been 
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Figure 7. (a: left) Mosaicked color image of Cas A with WIRC, with CO 
(red) , K-continuum (green) and P/3 (blue) . The range of surface brightnesses 
measured from the diffuse structures are 0.31 - 2.79xl0 -5 for CO, 0.23 - 
1.4xl0~ 5 for K-continuum, and 1.2 - 29.0xl0 -5 erg s 1 cm 2 sr 1 for P/3. 
CO-excess regions (in red) are marked as boxes (the bottom box marks the 
location of the Minkowski knot), (b: right) The spectral energy distribution of 
two representative CO emitting knots: knotl (green) and knot2 (red) which 
show excess emission at 2.294 /im. The contribution from the synchrotron 
continuum is shown by the dotted line. Comparison with the SED (scaled 
to the CO knotl continuum) of the NW ejecta knot (blue dashed curve) is 
shown. The SN 1987 spectrum is shown for comparison (orange curve; Mcikle 
et al. 1989). 



within a few years after the initial explosion. Our CO detection suggests that 
either the CO layer is not macroscopically mixed with ionised helium or that 
the helium is not ionised. Typical electron density of C- and O-rich regions dur 



ing f or mation of CO is show n to be 5xl0 8 cm 3 ( Gear hart. Wheeler fc Swart2 



1999). iTravaglio etaU (|l999l ) assumed that mixing occurs at the molecular level, 



prior to dust condensation, while Clayton & Nittler (2004), Deneault, Clayton 
&l Heger (2003) suggest that gases are unlikely to be mixed at the molecular 
level within a few years. Our CO detection suggests that molecular mixing is 
small (at least smaller than previously thought) in the SN ejecta and during the 
development of the reverse shock. 

7. Conclusion We present unambiguous evidence of dust formation in young 
SNRs. The dust map is almost identifical to the ejecta map, showing that dust 
is formed in the ejecta. Identified dust composition in SNRs include the features 
of Si0 2 [21 and 9 /im], MgSi0 3 [17 and 9.5 /mi], Si [15 /im], SiC [11.3 H and 
the featureless dust of carbon, AI2O3, and Fe. We also show that dust and ejecta 
compositions are closely correlated. E0102 shows dominant ejecta of Ne and O, 
while Cas A shows Ar, Si and Fe ejecta in additional to Ne and O. A total mass 
of freshly formed dust is 0.02- 0.054 M for Cas A, and 0.015 M Q for E0102. 
Gil. 2 also shows similar order of dust mass including Fe dust. The dust mass 
is more than one order higher than previous estimates but slightly lower than 
theoretical predictions. CO molecule detection from Cas A indicates on-going 
processes of astro-chemisty. 
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